Composite Fermions with Orbital Magnetization 
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For quantum Hall systems, in the limit of large magnetic field B (or equivalently small elec- 
tron band mass m,b), the static response of electrons to a spatially varying magnetic field is largely 
determined by kinetic energy considerations. This response is not correctly given in existing ap- 
proximations based on the Fermion Chern-Simons theory of the partially filled Landau level. We 
remedy this problem by attaching an orbital magnetization to each fermion to separate the current 
into magnetization and transport contributions, associated with the cyclotron and guiding center 
motions respectively. This leads to a Chern-Simons Fermi liquid description of the v = ^ state 
which correctly predicts the mt dependence of the static and dynamic response in the limit mt, — » 0. 



A useful theoretical tool for studying quantum Hall 
states is the composite fermions picture, known also as 
the Chern-Simons fermionic theory Based on the 
work of Jain Q , the Chern-Simons fermion picture was 
introduced by Lopez and Fradkin Q to study fractional 
quantized Hall states, and by Halperin, Lee and Read 
(HLR) as well as Kalmeyer and Zhang Q, to study 
even denominator filling fractions. 

In the composite fermion description, a system of spin 
polarized electrons in a two dimensional electron gas is 
transformed exactly into a system of composite, or gauge- 
transformed, fermions. These composite fermions carry 
an even number (j) — 2m of 'Chern-Simons' flux quanta 
and interact with each other both electrostatically and 
via the Chern-Simons gauge field 0. At a mean field 
level the fermions are subject to a uniform effective mag- 
netic field AB = B — &Qn e (f) where $o = h/e is the 
flux quantum, and n e is the average electron density 
(the speed of light c — 1). Such an approximation maps 
electrons at filling fraction v = jp- to fermions in zero 
magnetic field, where the ground state should be a gap- 
less Fermi liquid-like state. Similarly, this transforma- 
tion maps fractional quantized Hall states at filling 
to integer quantized Hall states of 



fractions v = „ p , 
fermions at v — p. 

The main systematic attempts for going beyond mean 
field theory have so far involved perturbative treat- 
ments of the Chern-Simons and electrostatic interac- 
tions. Attention has focused on one difficulty - 
the divergence of the composite fermion's effective mass 
at the Fermi surface when the range of electrostatic inter- 
action is shorter than, or equal to, that of the Coulomb 
interaction. This divergence, due to infra-red gauge field 
fluctuations, is manifested in the energy gaps of fractional 
quantized Hall states at v = 2r ^ p+1 (for large p) but does 

not affect the electronic linear response at v = Jj- at zero 
temperature, due to a mutual cancelation with another 
singular term Consequently @ the low energy ex- 

citations at v — 2p are best characterized by another, 
finite, effective mass, denoted by to*, which is the effec- 
tive mass of relevance to the present work. It is this to* 
which should determine the scale of the fractional Hall 



gaps for small values of p. 

In this paper we consider a second source of difficulty 
in the Chern-Simons theory which appears even when 
gauge-field fluctuations are not infra-red singular. (For 
example, there is no predicted infra-red divergences in 
the effective mass if the electron-electron repulsion falls 
off more slowly than 1/r.) This second difficulty is en- 
countered in the limit of a small band mass to;, or equiv- 
alently the large B limit. The fact that, in this limit, 
the electronic ground state and low energy excitations 
are constrained to the lowest Landau level, leads to cer- 
tain features of the electronic response to an external 
static vector potential which are not properly embod- 
ied in previous approximation schemes. In this paper we 
suggest a remedy to that flaw. We argue that while previ- 
ously used approximations are appropriate for calculat- 
ing the density-density response function Kqq (defined 
below) they are not appropriate for response functions 
such as the density response to a transverse vector po- 
tential Kqi. The approach we present allows us to derive 
an approximation to the linear response function which 
yields the correct to& dependence in the limit to& — > 0. 

Our approach is based on a separation of the current 
into a magnetization current which is associated with the 
cyclotron motion of electrons and a transport current as- 
sociated with the guiding center motion. This separation 
is achieved by attaching a magnetization /i M to each par- 
ticle. This magnetization originates from the electrons' 
orbital motion and is unrelated to the spin (which is as- 
sumed to be a frozen degree of freedom). In the limit 
TOb — > 0, the magnetization fi u is given by the Bohr 
magneton /if, = eS/(2mj). Our separation procedure, 
combined with approximations similar to those made in 
previous works, results in response functions that cor- 
rectly describe the to^ — > limit. 

We will mostly discuss the behavior of the electromag- 
netic response matrix which is closely related to the 
conductivity To define K, a weak vector potential 

A° xt is externally applied to a system at wavevector q 
and frequency uj, and consequently, a current is in- 
duced (Here Aq is the scalar potential, and jo is the in- 
duced density). We write the response function in the 
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form j fl (q,uj) = K \ iV (q 1 Ld)A c J Lt (q,Lo) where /i and v take 
the values 0, x, y. We will use the convention that the 
perturbation is applied with g||x so that the longitudi- 
nal current is j x = (u>/q)jo. Using the gauge A x = 0, 
we can then treat as a 2 x 2 matrix with indices 
taking the values or 1 denoting the time or transverse 
space components. In this notation the current vector 
iju is (io)i«)j and the vector potential is (Ao,A y ). In 
the HLR approach Q, the response matrix K is calcu- 
lated within the RPA. As reviewed below, an improved 
'Modified' RPA (MRPA) was introduced in later work 

To understand why we believe the RPA and MRPA 
miss some important physics, we consider a system at 
v = h such that the magnetic field is Bi/ 2 = 2n e <I>o. We 
now apply a weak static field 5B(r) at some small but 
finite wavevector q. At a mean field, RPA, or MRPA 
level, the transformed fermions see only the additional 
field 6B(r), and the energy is the same in regions where 
5B > as it is in regions where 5B < 0. However, the 
original electrons see the full field Bi/ 2 + SB(r). Since 
the electrons are in the lowest Landau level they have 
a kinetic energy ^hcu c (r) = fib(Bi/ 2 + SB(r)) which is 
lower in regions where SB < 0. In the limit of mj, — > 0, 
the electrons are then strongly attracted to regions where 
SB(r) < 0. This behavior should be reflected in Kqi, but 
is not accounted for at a mean field or (M)RPA level. 

In order to model this attraction to regions of minimal 
magnetic field, we attach to each particle a fixed magne- 
tization /x M which interacts with the magnetic field and 
creates an energy cost /j M B(r) as a function of position. 
In the limit where the interaction energy goes to zero 
(or equivalently if mi, — > 0), fi M is given by /if,, so that 
this energy cost consists of the kinetic energy only and 
is given by (i M B(r) — ^huj c (r). For finite interaction 
energy the local energy cost is composed of kinetic and 
interaction contributions. Within our framework, most 
of that energy should be described as magnetization en- 
ergy, with the remaining part described as quasiparticle 
interaction energy. There is, however, some arbitrariness 
in that separation, and consequently in the definition of 
/i M . For the present purpose we need to consider only 
the leading behavior in 1/mi, and we can therefore put 
/i M = fif,- Following the attachment of magnetization, 
we define an effective scalar potential seen by the mag- 
netized particles as Aq B = A + fi M B(r) and declare that 
the particles respond to this effective field. 

Another effect that is not accounted for in previous 
works involves the magnetization current, whose physi- 
cal interpretation is the following. In a high magnetic 
field, each particle can be thought of as traveling in a 
cyclotron orbit. When the density of particles is uni- 
form, the local currents of these orbits cancel and there 
is no net current in the bulk of the system. However, 
when there is a density inhomogeneity, these local cur- 
rents do not quite cancel and a net magnetization cur- 
rent results. Formally, we write jmag = V x M where 
M is the magnetization density which is given, for non- 



interacting particles in the lowest Landau level, by M = 
±zhu) c (r)n(r)/B(r) = fJ, b n(r) (where n(r) = n e + j (r) 
is the local density). Again, more generally /z& might be 
replaced by /i M such that M = /i M n(r). Thus, the mag- 
netization current is jmag = Mm(z x Vn(r)). By attaching 
magnetization fi m to each particle, which is equivalent to 
attaching a current loop, we separate out the contribu- 
tion of the magnetization current to the total current. 
It is then convenient to define the transport |J current 
jtrans = jtotai - jmag- The transport current may be inter- 
preted as the current of the bound particle-magnetization 
composite objects. 

When projected to the lowest Landau level, the pro- 
jected current and density operators satisfy PjP = 
Hb(Z x VPnP), ie, the projected current operator de- 
scribes only magnetization current ]To[ | (Here, P is the 
projection operator). Therefore, in the limit nib — > 0, if 
we apply a static scalar potential AQ Xt (q) to the system 
and we look at the current response in powers of l/m^, we 
find a transverse magnetization current fibZ x iqifoo^o 
originating from lowest Landau level matrix elements, as 
well as a current independent of m b , which may be de- 
scribed as arising from virtual transitions to other Lan- 
dau levels. Thus, if q is finite and mb — > 0, we expect 
Kiq/Kqq — > iq/ib- This result is not contained in previous 
works based on the Chern-Simons approach. 

Before detailing our remedies to the above problems, 
we review the RPA and MRPA. Both approaches sep- 
arate the long ranged interactions by writing K as 

K-i=TT* + U where U = ( ^ ~f> ) (1) 

Here v(q) — e 2 /(eq) is the Coulomb interaction and 
c(q) = icp&a/q is the Chern-Simons interaction. In other 
words, n is the part of K that is irreducible with respect 
to both v(q) and c(q). The RPA is obtained by 
setting n equal to K°, the response of a noninteracting 
system of fermions with mass mb in the mean magnetic 
field AB. This RPA approach has the problem that the 
low energy excitations are on a scale set by the cyclotron 
energy (ie, by mb) rather than by the interaction energy. 
The MRPA || repairs this problem by defining an ef- 
fective mass m* that is set phenomenologically by the 
interaction scale 0. Galilean invariance is maintained 
by adding a Landau interaction term T\ to insure that 
Kohn's theorem and the /-sum rule will be satisfied. To 
define the MRPA we write 

n- 1 = [n*]- 1 + Tx (2) 
(m*-m & ) /^ o \ 

n e e 2 \ -1 J y ' 

The MRPA is then obtained by setting n* equal to the 
response K°* of a system of noninteracting fermions of 
mass m* in the mean magnetic field AB. The response 
function thus calculated (using n* = K°* and Eqns. || 
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and P will be called K MB - FA . Comparisons of numerical 
results of exact diagonalizations to results of i^oo calcu- 
lated in the MRPA were quite favorable [|| for the low 
energy excitations at v = 2m P p+i ^ or sma H P- 

As described above, the (M)RPA approach does not 
account for magnetization effects due to the fact that 
when we take the mean field solution as a starting point 
for a perturbation theory for the Chern-Simons fermions, 
we lose the fact that the original electrons travel in local 
cyclotron orbits. This piece of physics would presumably 
be recovered if one could correctly carry out a renormal- 
ization procedure which eliminates the high frequency 
motion at the cyclotron frequency to obtain equations 
valid on the scale of the Coulomb interaction. In the 
current work, we recover this physics by attaching mag- 
netization to each particle by hand. This attachment 
is not an exact transformation, but is rather a way of 
modeling behavior that is lost when we take the mean 
field as a starting point. As we will see below, within 
a Landau-Fermi liquid theory picture, this attachment 
seems to give the correct quasiparticles for the system. 

The attachment of magnetization to the particles does 
not change the interparticle interactions. However, the 
magnetized fermions now respond to the effective poten- 
tial and the motion of these magnetized fermions yields 
only the transport current response. We thus define a 
matrix K to be the transport current response of the elec- 
trons to the external effective potential. We thus write 



A cfi = MU 

jtotal = -Mjtrans 



M 



1 

iqu, M 1 



(4) 



and relate K to K by K = MKM^ . As discussed above, 
in the limit nib — > 0, we must have ji M — > /it,. In the rest 
of this paper, however, we will consider ji M = 

We now propose an approximation which we call the 
'Magnetized Modified RPA' or M 2 RPA where we set K 
equal to X MRPA . Thus in the M 2 RPA we approximate 



(5) 



This equation is a central result of this paper. In the 
limit nib — * 0, the M 2 RPA correctly describes the static 
response properties described above. It should be noted 
that the M 2 RPA predicts K q = K^ a RFA and is there- 
fore equally well supported by exact diagonalizations || . 
Furthermore, we note that the excitation modes of a sys- 
tem are determined by the poles of det[JC], and since 
det[K] = det[K], the MRPA and the M 2 RPA predict the 
same excitation energies. However, the M 2 RPA calcula- 
tions of K i and K w differ from K^ RPA and Kf Q RPA by 
terms of order q/mb and the M 2 RPA value of K\\ dif- 
fers from Ki 1RFA by terms of order q/nib and by terms 
of order (q/nib) 2 - It should be noted that all finite q ex- 
perimental tests H of the Chern-Simons theory to date 
have measured only A'oo- Similar to MRPA, we expect 
the M 2 RPA, in addition to describing the v = ^ Fermi 
liquid states, should properly describe the Jain series of 



quantized states v = 2r Jp+i f° r sma U P- At large values 
of p, the description should be modified to account for 
the effects of the singular infra-red gauge fluctuations. 

We now turn to discuss how the M 2 RPA fits into the 
general picture of a Fermi liquid theory of the v = l/(2m) 
state. In essence, we show that M 2 RPA amounts to 
adopting the Fermi liquid picture of Ref. M as describing 
the dynamics of magnetized composite fermion quasipar- 
ticles rather than unmagnetized ones. 

In Fermi liquid theory for fermions with short ranged 
interactions, such as 3 He, the response function K is de- 
scribed by the solution of a Landau-Boltzmann equation 
lyj . However, for fermions with long ranged interac- 
tions jllj , the Silin extension of the Landau theory as- 
serts that it is not the function K that is described by 
the Landau-Boltzmann equation, but rather a response 
function called n. In the case of electrons at zero mag- 
netic field in 2D, for example, n is defined by Eq. [j], but 
with c(q) = 0. Eq. [I] separates out the Hartree part of 
the diverging long ranged interaction such that n gives 
the quasiparticle response to the sum of the external vec- 
tor potential and the induced internal vector potential. 
Within this theory, the Landau-Boltzmann equation de- 
scribes the dynamics of quasiparticles of mass m* near 
the Fermi surface The short ranged effective 

interactions between two quasiparticles whose momenta 
differ by an angle 8 is described by a function f(9). The 
'Landau coefficients' are defined as // = ^ J Q * d0f{9)e lW 
(by symmetry, /; = /_;). The parameters fa and /i are 
fixed by the identities jll| 
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+ /o 
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nib 



h 



m ■ 



2nh 2 . 



(0) 



The second identity is a result of Galilean invariance [^l) . 

Since in the Chern-Simons problem, the interactions 
are also long ranged, we have similarly made the separa- 
tion defined in Eq. [I] (with c(q) — i$o(j)/q). However, 
in the limit nib -> 0, further separation should be carried 
out to remove the diverging magnetization effects. To 
this end we define a response function n by 



n = MUMl 



(7) 



definition, n relates the transport current of the magne 
tized quasiparticles to the effective vector potential (in- 
eluding both external and internally induced contribu 
tions ]l0[ ] ) . It is A which we claim is given by a Landau- 
Boltzmann equation describing the dynamics of quasi- 
particles with the finite effective mass m*. 

Recent works attempting to formulate a Landau- 
Boltzmann equation for the Chern-Simons problem 



pointed out that in Eq. ||, ^ is taken at fixed AB. This 
yields an /o (as well as an fx) on the scale of ©(m^ -1 ) 
in the limit nib — * 0. We will now show how an O^m^ 1 ) 
value of fo is consistent with our attachment of magne- 
tization. Fixing AB at zero means that when we add 
a fermion to the system, the external field must also be 
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increased by a total of <fi flux quanta. Thus, the magnetic 
field is linked to the density n via B = 0n$o- The inter- 
action energy between the magnetization M = U(,n and 
the external field is given by E — M • B = TT(f)Ti 2 n 2 /nib. 
Differentiating this with respect to n we obtain a mag- 
netization contribution to the chemical potential /i mag = 
2tt (f>h 2 n/mb — fiuj c such that the magnetization contri- 
bution to fo is fo = d/z mag /dn = 2n(f)h 2 /nib which is also 
the inverse compressibility of free electrons of mass nib a t 
constant AB. The coefficient f is written f = fo + Sfo 
where fo is Oijn^ 1 ) and Sfo is on the smaller interaction 
scale. As mentioned in Ref. 0, in the limit nib — ► 0, 
the requirement that the low energy spectrum is inde- 
pendent of TOfc forces the other interaction coefficients (// 
for / ^ 0, 1) to be on the interaction scale Jl0[ . 

Since in the limit of m& — ► 0, fo and f\ are on the 
bare mass scale whereas all other coefficients /; (as well 
as Sfo) are on the interaction scale, we will separate out 
the contributions of these two coefficients by writing |1G] 



rr 1 = [riT 1 + f + Ti 



(8) 



where To = diag[/o,0]. The function II* is to be calcu- 
lated using a Landau-Boltzmann equation representing 
quasiparticles with the same effective mass m* and in- 
teraction coefficients fi except that f\ is artificially set 
to zero and the magnetic contribution fo is subtracted 
off of fo- The separation of the coefficient fo, equivalent 
to taking v(q) — ► v(q) + fo in Eq. [|, is achieved by not- 
ing that fo corresponds to a short ranged density-density 
interaction. The separation of the nonzero f\ coefficient 
P,p"o| is analogous to that described in Eq. |2| and is de- 
rived in Ref. j6| . Having made this separation, we expect 
that the response H*(q,u>) is independent of nib in the 
limit nib —> and is well behaved for all values of q/nib- 

In Fermi liquid theory, the Landau-Boltzmann equa- 
tion does not correctly describe the Landau diamagnetic 
contribution to the transverse static response. Similarly, 
we suspect that here the function 11*^ derived from the 
Landau-Boltzmann equation lacks a term of the form q 2 \ 
where x is some appropriate Landau susceptibility which 
we expect to be on the scale of the interaction strength. 
As usual, if we fix the ratio to / q to be nonzero, and take 
q — > 0, this diamagnetic term becomes negligible. 

Using the identities To = M^TqM and U + To = 
h'^^UM- 1 (which holds in the limit m b -> 0), we find 
that M 2 RPA is equivalent to approximating II* by K°*, 
the response of a free Fermi gas of particles of mass m* , 
and calculating the response using Eqns. |l|, fj], and ||. 

While the M 2 RPA describes well the excitation spec- 



trum of fractionally quantized Hall states v 



for 



2mp+l 

large p at low wavevector q, it does not do so at high q. In 
particular, excitations at high q are sensitive to the infra- 
red divergence of the effective mass due to the gauge field 
fluctuations which are neg lected in 1VTRPA. 

To conclude, the transformation Eqns. [j], ||, ||, 0, and 
|8| do not in themselves involve any approximations, and 
may be considered simply as a means of defining a new 



'irreducible' response function Tl*(q,Lu). Our claim that 
II* is well behaved in the limit mj, — » we believe to 
be an exact statement (although we have not proved it 
rigorously) independent of the approximation used to de- 
fine the M 2 RPA. (For this to be true it does not matter 
whether we take the moment u M in Eq. ^ to be precisely 
fib or whether we include a correction of order fibnib/m* .) 

The approximation we introduce, the M 2 RPA, de- 
scribes the v = state as a Fermi liquid of magnetized 
composite fermions with a finite renormalized effective 
mass m* , an f\ parameter dictated by Galilean invari- 
ance and an fo parameter originating from the interac- 
tion of the magnetization with the magnetic field. All re- 
maining Fermi liquid parameters (which are expected to 
be on the much smaller interaction scale) are neglected. 
The M 2 RPA predicts the same Kqo as the MRPA, but in 
contrast it yields the correct behavior for Koi,Kio, and 
Kn in the limit nib — * for arbitrary small q. In the 
limit q -> 0, for fixed m b ^ 0, the M 2 RPA and MRPA 
become identical for all components of K^ v . 
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